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Abstract 

Aretakis has proved that a massless scalar field has an instability at the horizon 
of an extreme Reissner-Nordstrom black hole. We show that a similar instability oc- 
curs also for a massive scalar field and for coupled linearized gravitational and elec- 
tromagnetic perturbations. We present numerical results for the late time behaviour 
of massless and massive scalar fields in the extreme RN background and show that 
instabilities are present for initial perturbations supported outside the horizon, e.g. an 
ingoing wavepacket. For a massless scalar we show that the numerical results for the late 
time behaviour are reproduced by an analytic calculation in the near-horizon geometry. 
We relate Aretakis' conserved quantities at the future horizon to the Newman-Penrose 
conserved quantities at future null infinity. 
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1 Introduction 



Supersymmetric black holes play an important role in string theory. Given their impor- 
tance, it is natural to ask whether or not they are classically stable: does a small initial 
perturbation remain small under time evolution? 

In a supergravity theory, the fact that a supersymmetric solution saturates a BPS bound, 
and therefore minimises the energy (at fixed charge), does not imply classical stability. For 
example, the classical stability of Minkowski spacetime in vacuum GR does not follow from 
the positive energy theorem. Instead, the proof involves a lengthy analysis of the Einstein 
equation |T]. Furthermore, Anti-de Sitter spacetime is a supersymmetric solution of various 
supergravity theories but nevertheless it is classically unstable against the formation of small 
black holes [2]. 

In these examples we are referring to classical stability under time evolution determined by 
the nonlinear Einstein equation. Recently, Aretakis has proved that even linear perturbations 
of a supersymmetric black hole can exhibit an instability [3l He considered an extreme 
Reissner-Nordstrom black hole. In ingoing Eddington Finkelstein (EF) coordinates, the metric 



The future event horizon Ti'^ is at r = M. The generator of time translations is d/dv. This 
solution preserves half of the supersymmetry in minimal N = 2 supergravity [6]. 

Aretakis studied a massless scalar field ip in this spacetime. Consider initial data for ip 
specified on a spacelike surface S intersecting "H^ and extending to infinity (Fig. [1]), with 
i(j decaying at infinity. This uniquely determines ip in the future domain of dependence of 
S, which contains the part of the black hole exterior that lies to the future of S, and a 
neighbourhood of the part of 'H'^ that lies to the future of S. Aretakis proved that the 
solution exhibits both stable and unstable features. Stable features are that ip decays at late 
time (large v) on, and outside, and all derivatives of decay outside T-L^. However, by 
constructing certain conserved quantities on "H"*" (reviewed below), Aretakis proved that drijj 
generically does not decay at late time on l-C^ . This implies the existence of an instability: 
since dr'ip decays for r > M but not for r = M, it follows that d^ip must blow up at late time 
on Aretakis proved that d^tp generically blows up at least as fast as v''~^ at late time on 



This instability is not a coordinate effect because —d/dr can be invariantly defined as the 
tangent to ingoing radial null geodesies of unit energy. The instability involves polynomial 
growth in time, which would make it hard to discover using a mode analysis. The above 
results imply that the component T^r of the energy momentum tensor of ip decays at late 
time outside "H^ but not on l-i^ . One might regard this as "hair" on the horizon of the black 
hole. Note that T^^ is closely related to the energy density measured by an infalling observer. 

We emphasise that non-extreme black holes such as Schwarzschild or non-extreme Kerr 
have been proved to be stable against massless scalar field perturbations: ip and all its deriva- 
tives decay on, and outside, T-C^ [71 18]. One might think that extreme black holes also should 

^ Ref. [5] conjectured the existence of an instability of extreme black holes. 
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Figure 1: Penrose diagram for extreme RN black holes. Aretakis took initial data specified 
on a spacelike surface S intersecting "H"*" and extending to infinity. 

be stable because a generic perturbation would make the black hole non-extreme. This is 
incorrect. First, this involves backreaction on the metric, which is a nonlinear effect. If the 
scalar field perturbation is 0{e) then the resulting metric perturbation will be O(e^). A sec- 
ond order effect won't eliminate a linear instability. What might happen is that the endpoint 
of the instability is generically a non-extreme black hole, i.e., derivatives of ip become large 
along "H"*" (so there is still an instability) but eventually are damped by nonlinear effects. 
Second, even though generic initial perturbations will make the final black hole non-extreme, 
there are probably non-generic initial perturbations for which the horizon is extreme at late 
time. 

Minimal N = 2 supergravity does not contain a scalar field but the Aretakis instability 
can be embedded in a supersymmetric theory as follows. Consider a supersymmetric static 4- 
charge black hole solution of type II supergravity compactified on T^. In general, this solution 
has non-trivial moduli fields but if the 4 charges are set equal then the moduli are all constant 
and the geometry is that of extreme RN. Linearized fluctuations in the moduli are massless 
scalars, hence Aretakis' result implies the existence of a linearized instability at the horizon 
of an extreme RN black hole in this theory. 

Aretakis has considered also the case of a massless scalar field ip in the extreme Kerr 
geometry. He has proved decay of an axisymmetric field [9] and that derivatives of ip exhibit 
an instability at H"*" analogous to the extreme RN instability ^U\. Ref. [H] generalised these 
results to show that a massless scalar field instability occurs at for any extreme black hole 
and that a similar instability occurs for linearized gravitational perturbations of the extreme 
Kerr solution. Ref. [12] has extended the latter result to linearized gravitational perturbations 
of certain higher-dimensional extreme vacuum black holes. 

Aretakis' proof of instability (reviewed in section 12. 2p involves an infinite set of conserved 
quantities at 'H'^, linear in ip. We will call these the Aretakis constants. We will show in section 
12.31 that these constants are closely related to a set of conserved quantities at X"*" (future null 
infinity): the Newman- Penrose constants [13]. Indeed, there is a conformal isometry of the 



4 



extreme RN geometry which interchanges T-C^ with X"*" [Hj. This map exchanges the Aretakis 
constants with the NP constants. 

The Aretakis instabihty is associated to outgoing radiation at "H"*": the proof of instabihty 
requires that the initial data be non-vanishing at and, as we will explain below (section 
12.11) . the instability is closely related to the absence of a redshift for outgoing photons at "H"*". 
For string theory applications, it seems more natural to consider what happens if one perturbs 
an extreme RN black hole by dropping something into it. Can one trigger an instability using 
ingoing radiation? Some stability results for this case were presented in Ref. [T3], which 
considered initial data for iIj compactly supported outside l-i^ . It was shown that the scalar 
field and all of its derivatives remain bounded in static coordinates. However, such coordinates 
cover only the black hole exterior so this result does not exclude the existence of an instability 
at H+. 

We have investigated this problem numerically and describe our results in section 121 We 
find that initial data corresponding to an ingoing wavepacket does lead to an instability at 

. It afflicts quantities with one more r-derivative than in the outgoing case: ip and dril^ 
decay on, and outside "H"*", d^ilj decays outside "H"*" but generically does not decay on "H"*" 
and dlip generically blows up at late time on "H^. We also study the late time behaviour 
(tail) of ip and find that the field decays more slowly if the Aretakis constants are non-zero. 
This extends previous numerical work on massless scalar field tails in extreme RN [TB], which 
considered outgoing wavepackets outside 'H'^, which have vanishing Aretakis constants. 

It is natural to ask whether the Aretakis instability can be seen in the AdS2 x near- 
horizon geometry of extreme RN. In section H] we will show that it can: it occurs at the horizon 
of AdS2 in Poincare coordinates. This is not in contradiction with known stability results for 
linear fields in AdS2 because it turns out to be a coordinate effect: there is no invariant way 
of defining coordinates in AdS2 analogous to the (f,r) of ([1]). Nevertheless, the AdS2 results 
are interesting because they are in excellent agreement with our numerical results for the late 
time behaviour of ip at 'H'^ for the extreme RN spacetime. 

The remainder of our paper concerns generalisations of Aretakis' work to other fields. 
In section \5\ we consider a massive scalar field. For discrete values of the mass m (m^ = 
n{n + 1)M~^ for positive integer n) we prove the existence of conserved quantities analogous 
to the Aretakis constants and use these to prove instability at T-L~^. For more general values 
of the mass we demonstrate instability numerically. A more massive field is more stable: the 
number of r-derivatives exhibiting decay at "H"*" increases with the mass of the field. However, 
our numerical results indicate that for any mass there exists k such that d^ip generically blows 
up at late time on "H"*". 

Finally, it is desirable to have an argument for instability of extreme RN that does not 
rely on the existence of a scalar field. In section [6] we will consider (coupled) gravitational and 
electromagnetic perturbations. We prove that for all types of perturbations of this kind, there 
exist towers of conserved quantities on the horizon, which can be used to exhibit instabilities 
very similar to those of the massless scalar field. We construct an explicit gauge invariant 
combination of the Maxwell field strength perturbation, the metric perturbation and its first 
derivatives. The instability is strongest for the I = 2 multipole moment, for which we show 
that a certain combination of the first and second r-derivatives of this quantity is conserved 
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on and hence generically cannot decay. We argue that this imphes that a quantity with 
one more r-derivative generically blows up at late time on "H"*". Hence there is an instability 
at the horizon of extreme RN in Einstein-Maxwell theory (or minimal N = 2 supergravity). 

Note added: As this paper was nearing completion, Refs. [211 132] appeared. Ref. |3I] has 
significant overlap with our section 12.31 relating the Aretakis and NP constants. Motivated 
by our numerical results, Ref. [32] gives a proof of the existence of an instability for ingoing 
massless scalar field radiation. Ref. [31] also gives an argument for the existence of such an 
instability. 

2 Massless scalar horizon instability 

2.1 Motivation for instability 

To understand why there might be an instability at the horizon of an extreme black hole, 
it is useful to consider first the case of a massless scalar field if) in a, non-extreme black hole 
spacetime. A first step in trying to prove stability is to consider the energy of ip, written as an 
integral -EfS] over a spacelike surface S extending from the future horizon to infinity. This 
is a non- increasing function of time, i.e., < if E' lies to the future of E. Outside 

"H"*^, the integrand, i.e., the energy density, is a positive definite function of d^j^ip. Hence the 
fact that E is non-increasing implies that d^ip cannot become large. However, precisely on 
"H"^, the energy density degenerates: it does not depend on the derivative of if) transverse to 
T-i^ . Hence it is consistent with energy conservation for this derivative to behave badly. 

For non-extreme black holes, such behaviour has been excluded using the horizon redshift 
effect. This is the fact that the energy of a photon moving tangential to undergoes a 
redshift proportional to e~'^^ where k, is the surface gravity and v the Killing time along the 
horizon. Using the wave analogue of this effect is an important step in the proof that 
and all of its derivative decay on and outside Ti'^ in a Schwarzschild [7| or non-extreme Kerr 
spacetime [8]. 

Now consider an extreme black hole. Such a solution has k = 0: the horizon redshift effect 
is absent and so the problem of controlling derivatives transverse to 'H'^ remains. Aretakis 
has proved that this problem cannot be overcome for extreme RN and extreme Kerr. Ref. 
[TT] showed that it cannot be overcome for any extreme black hole. Below we will review 
Aretakis' argument for extreme RN. 

We note that Ref. [17] considered the case of a massless scalar field in RN and argued 
that there is a qualitative difference between the extreme and non-extreme cases arising from 
the behaviour of outgoing waves near T-L~^. However, the detailed predictions of Ref. [T7| are 
in disagreement with numerical results of Ref. [TB] and also with the numerical results that 
we will present below. 
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2.2 Aretakis' argument 

The equation of motion of a massless scalar is 

V^^/^ = . (2) 
Working in the coordinates of ([1]), first expand ip in spherical harmonics: 

oo 

i;{v,r,n) = Y,U'"^r)Yi{n) (3) 

(=0 

(we suppress the azimuthal index m) and substitute into the equation of motion to obtain 

2rdMr^lji) + dr {^drilJi) - l{l + 1)^^/ = (4) 
where A = (r — M)^. Set / = and evaluating at "H^ (r = M) shows that 

i7oM = ^[9r.(rV^o)].=M (5) 

is conserved, i.e., independent of v. 

For generic initial data, Hq will be non-zero. Hence it remains non-zero. Therefore ijj and 
drip cannot both decay at 'H'^. In fact, ip decays at late time [3] hence it follows that its 
transverse derivative at 'H'^ does not decay: 

{dr'lpo)r=M Hq as V OO . (6) 

The rr component of the energy momentum tensor of ijj is Trr = {drip^, which generically 
does not decay at late time on 'H'^. This implies that the energy density measured by an 
ingoing observer at T-L^ does not decay, as suggested by the absence of the horizon redshift 
effect. 

Now act on (jl]) with dr, set / = and evaluate at "H"*" to obtain 

[d,d',{riPo) + driPo]^^^j = 0. (7) 

Hence 

[d^d^{rilJo)]^^j^j ^ -Ho as v ^ oo (8) 

which, together with decay of ipo and (E]), implies blow-up of the second transverse derivative 
of ip at late time on "H^: 

{d^ipo)r=M -j^v as V ^ oo . (9) 

By taking further r-derivatives of (j4]) it can be shown that {d^ipo)r=M oc w^"^ for large v [4]. 

This instability is a property of s-wave (/ = 0) perturbations. There is a corresponding 
result for / > 0: acting on (jl]) with dl and evaluating at r = M reveals that 

Him^^HlrdririPi)]}^^^ (10) 

is conserved. Aretakis shows that d'^ipi decays on, and outside 'H'^ for /c < / [3]. Using this, it 
follows that d^y^ipi generically does not decay at and, arguing as above, d^'^^ipi generically 
blows up linearly at "H"*" 
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2.3 Conformal isometry 

The extreme RN solution has a discrete conformal isometry [H]. In static coordinates 
(t, r, 9, 0) it is 

$ : (t, r, e, 0) ^ (t, r' = M + -, 6, 0) (11) 

r — M 

This is self-inverse: $ = The push-forward of the metric under this diffeomorphism is 

Mg) = ^'g (12) 

where 

M 

Q = — — . (13) 
r-M ^ ' 

To understand how this acts on the horizon, we will write it in Eddington-Finkelstein coor- 
dinates, related to the static coordinates as follows. First define the tortoise coordinates 

r-M\\ 



r,(r) = r-M + 2M\n [^^j^ j " , (14) 

and note that $ sends r* to r[ = — r*. Now let 

u = t - r^, V = t + r^ (15) 

and so $ sends u to u' = v and v to v' = u. Hence $ maps the point with ingoing EF 
coordinates {v,r,6,(j)) to the point with outgoing EF coordinates {u' = v,r',6,(f)). It follows 
that $ maps 'H'^ to and vice-versa. 

The extreme RN geometry has vanishing Ricci scalar, which implies that the massless 
scalar wave equation is conformally covariant: if is a solution in the extreme RN geometry 
with metric g then VL'^ip is a solution in the geometry with metric Vt^g. Hence the conformal 
isometry can be used to generate a new solution if) of the massless scalar wave equation from 
an old one = Qtp o $. In coordinates: 

M 

^{u, r, e, 0)o = ;—]^Hv' = r', 9, <P)i (16) 

On the left, the subscript O indicates that outgoing EF coordinates are used, on the right the 
subscript / indicates that ingoing EF coordinates are used. 

What happens to the Aretakis conserved quantities under this map? Long ago Newman 
and Penrose (NP) argued that there is an infinite set of conserved quantities associated to 
linear massless fields at X"*" [13]. These can be defined by developing an asymptotic expansion 
in inverse powers of r near X"*" : for the / = multipole moment of ip assume the expansion 

Wo ~ \ 5— (1 ' J 

as r — 7- 00 1^ The equation of motion for il) implies that the quantity /i is independent of 
u: this is the NP constant of which we will denote by /if^/"]. It is easy to check that the 



^This can be deduced if we assume that VI ^tp is smooth at 1+ in the conformahy compactificd spacetime. 
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conformal isometry exchanges the Aretakis conserved quantity with the NP constant, i.e., the 
Aretakis conserved quantity for ipQ is the NP constant of ipQ and vice versa. Exphcitly we find 



Horn 



AM 



(18) 



M3 



The physical interpretation of the NP constants is not well-understood but it is clear that 
they are closely related to ingoing radiation at just as the Aretakis constants are related 
to outgoing radiation at H^. The NP constants influence the decay rate (tail) of outside 
the horizon: there is evidence that initial data with non-vanishing NP constants results in 
slower decay than data with vanishing NP constants |18]. We will return to this point below. 

As explained above, if any of the Aretakis conserved quantities of is non-zero then 
ip has an instability at T-L^. Does this mean that there is an instability of ip at X+ if its 
NP constants are non-zero? Even in Minkowski spacetime, there is a late-time blow-up of 
transverse derivative oiQ~^ip at X"*" when the NP constants are non-zero [13]. However, unlike 
"H"^, X"*" is not part of the physical spacetime hence this is not an instability: the field and all 
its derivatives decay in the physical spacetime. 

The definition of the NP constants has been criticized because of the assumed smoothness 
of at X"*" [in]. For spacetimes with non- vanishing NP constants, generically one expects 

only a finite degree of differentiability at X+. In contrast, ip will be smooth at 'H'^ if it arises 
from smooth initial data prescribed on a surface E which extends a finite distance behind 
as in Fig. [1] Thus there is an asymmetry present in the amount of differentiability to be 
expected at and X+. Starting from initial data specified on S we get a solution ip smooth 
at H"*" but generically not smooth at X+. Applying the above map then gives a solution ip 
for which Q~^ip is smooth at X+ but ip is not smooth at H"^. Such tp does not correspond to 
smooth initial data on S and therefore lies outside the class of solutions considered in this 
paper. 

3 Numerical results for massless scalar 

3.1 Double null coordinates 

A crucial assumption in Aretakis' proof of instability is that the conserved quantities Hi[ip] 
are non-vanishing. This requires that the scalar field ip be non-vanishing at the intersection 
of "H^ and the hypersurface S on which initial data is specified. Hence the results in section 
12.21 do not reveal what happens if one perturbs an extreme RN black hole by sending waves 
into the black hole from outside. We will investigate this problem numerically. To do this we 
first introduce coordinates which are better suited to numerical evolution. 

In [u, t')-coordinates the extreme RN black hole metric is 



ds 



2 



F{r{u,v))dudv + r{u,vYdfl'^ , 



(19) 



where 




(20) 
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and r{u,v) can be determined by solving r^{r) = (v — u)/2 with defined by ^ 

This metric is singular at . Since we wish to investigate time evolution of perturbations 
on the horizon, we need a regular metric there. Thus, we introduce a further coordinate 
transformation defined by 

^ = -r.(M-f/) = r/-2Ml„(M)-^. (21) 

From the definition of r*, we have du/dU = 2/F{M — U). In (t/, w)-coordinates "H^ is located 
at = and and f/ < corresponds to the exterior of the black hole. In this coordinate 
system the metric reads 

ds^ = '^£!f)-^dUdv + r^rffi^ , (22) 

F(M -U) 

where r = r{U,v). Using r*(r) = [v — u)/2 it can be seen that r is analytic in U,v. In 
particular we can expand r for small U (and fixed v) to obtain 

It follows that F{r)/F{M - U) = I + 0{U) for small U and the metric ([22]) is analytic at 
1 i.e., it can be analytically continued to the black hole interior [/ > 0. 

3.2 Klein-Gordon equation and initial data 

Consider the Klein-Gordon equation 

(V^ - m^)i) = (24) 

in the extreme RN background. Defining (p = rtpi, where ipi is the Ith multipole moment of 
ip as in equation (E]), we obtain a wave equation for 0. In the (f/, f )-coordinates introduced 
above this reads 

-Adud,<j) = V{U,v)<j) (25) 
,y,rr ^ 2F(r) fF'{r) /(/ + !) ^ 

In this section, we study the massless scalar field m = 0. We will consider the massive scalar 
in section [5l 

We consider a null "initial" surface defined by 

So = {f/ = Uo, v>vo}U{U> Uo, V = vo} (26) 

and impose the following initial data: 

mvo)=exp > 0(?7o,t;) = O, (27) 

10 




Figure 2: Schematic plot of the setup used in the numerical calculations. 

or ^ ^ 

m vo) = 0, 0(f/o, v) = exp (- ^""'jP ) , (28) 

which correspond to an outgoing and ingoing wavepacket, respectively. In Fig. [2] we give a 
schematic plot of the setup. 

We will solve this problem numerically. In our numerical calculations, we use units such 
that M = 1 and set Uq = —0.5 and Vq = 0. We describe the details of this in Appendix lA.ll 

3.3 Numerical simulation for spherically symmetric mode 

In this section we study spherically symmetric solutions, i.e., we set / = 0. In this case 
Aretakis' conserved quantity defined in Eq. (|5]) is simply given by Ho[ip] = dr4>\r=M- Hence 
the outgoing wave initial data (!27|) has Hq[iIj] ^ (unless /i = 0) whereas the ingoing wave 
initial data (l28l) has Ho[ip] = (since its support does not intersect the horizon). 

3.3.1 Non-zero Aretakis constant 

Firstly, we consider solutions for which -f/^of"^] 7^ 0. Although, in this case, the instability 
of the scalar field has been shown analytically, as reviewed in section \2.2\ we will investigate 
finer details of the time evolution of the scalar field by solving the wave equation numerically. 

Although we used the {U, v) coordinates for our numerical calculations, we will display our 
results using [v, r) coordinates. The reason for this is that the [v, r) coordinates are a preferred 
set of coordinates associated to the symmetries of the background. In particular, d/dv is the 
generator of time translations and —d/dr is tangent to ingoing radial null geodesies of unit 
energy. The {U, v) coordinates are not so closely associated to geometric invariants e.g. the 
generator of time translations is not d/dv in these coordinates. Hence when working in {U, v) 
coordinates it is less clear whether something is a physical effect or merely a coordinate effect. 
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Figure 3: Functions </)(y, r) and dr4>{v, r) for / = and i^o 7^ on fixed v slices. We consider 
outgoing wave initial data fl27j) with (ex, /.i) = (0.1, —0.1). The horizon is at r = 1. We can see 
that dr4> becomes steeper near the horizon as v increases. This implies that blows up on 
the horizon at large v. 
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V V 



Figure 4: Time dependence of 4>\r=M and d^(j)\r=M for / = and Hq ^ 0. We choose param- 
eters as (cr, /i) = (0.1, —0.1), (0.05, —0.1), (0.1, —0.05) in the outgoing wave initial data ( 1271) . 
We can see that d'^.(^\r=M blows up as ~ v. 

Consider outgoing initial data (!27|) with (a,/!) = (0.1, —0.1). The time evolution of and 
dr4> in (w, r)-coordinate is plotted in Fig. |31 We can see that decays as v increases. On 
the other hand, 9^0 does not decay along the horizon because dr(^\y=M must be conserved. 
However, outside the horizon, 9^0 decays. As a result, 9^0 becomes steeper near the horizon 
as time increases. This is consistent with the fact that 9^0|r=M must blow up along the 
horizon. 

In Fig. m we plot the time evolution of 0|r=Af and d1<^\r=M for various initial data: 
(a, yu) = (0.1, —0.1), (0.05, —0.1), (0.1, —0.05). We can see that the time behaviour of 9^0|r=Af 
is consistent with the expected linear blow up shown in section 12.21 In Appendix IA.2t we 
describe our method for evaluating the transverse derivatives on the horizon c?"0|r-=M- 

Now, we investigate the late time behaviour of 0|r=M- By fitting the absolute value of 
0|r=M to a power law decay u'^, we obtain the following exponents: a = —0.999, —1.018, —0.991 
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for (cT,/i) = (0.1, -0.1), (0.05, -0.1), (0.1, -0.05), respectively^ We used the numerical data 
in the range 1000 < v < 2000 for the fitting. These results suggest that the scalar field on the 
horizon (j)\r=M decays as at late time. To determine the coefficient of the power law decay, 
we fit the 4'\r=M to the function aHo/v + ^ for 1000 < v < 2000. Then, we find the fitting 
parameter a = -1.997 for all initial data: (a,/i) = (0.1, -0.1), (0.05, -0.1), (0.1, -0.05). 
Hence, for all cases that we have considered, the late time behaviour of 4'\r=M is given by 

2.00i/o . ^ 

(p\r=M ~ V OO . (29) 

V 

In Sec. Hlwe will show that the same expression can be obtained analytically from an AdS2 
calculation. In Ref. |4] it was shown analytically that the scalar field on the horizon at late 
times decays at least as fast as v~'^^^, which is consistent with our results. 



3.3.2 Zero Aretakis constant 

Now consider the case where Hq[xIj] = 0. This includes the case of an ingoing wavepacket 
as well as an outgoing wavepacket with n = 0. The proof of the Aretakis instability in section 
12.21 does not apply now; instead we will use a numerical calculation to study this case. For 
the outgoing wave (1271) . we take /i = and a = 0.05, 0.1, 0.15. For the ingoing wave ( l28l) . we 
choose the parameters {cr',^') = (3.0,10.0). In Fig. |5]we plot (p\r=M, d^(f)\r=M and 9^0|r=M 
as functions of v. Our results are consistent with the following behaviour: (j)\r=M decays, 
d^(f)\r=M approaches a non-zero constant and df(j)\r=M blows up proportionally to v. Thus, 
we conclude that there is an instability even for initial data with Ho[ip] = 0. 

Now consider the late time behaviour of 4>\r=M- Fig. [5] suggests that 4>\r=M decays as a 
power law v"" at late time. Fitting to such a decay law we find the exponents a = —2.000, 
—1.997, —1.995 for the outgoing wave with a = 0.05, 0.1, 0.15 and a = —1.975 for the ingoing 
wave. This suggests that, for Hq = 0, the late time behaviour of 4>\r=M is 0|r=Af ~ C/v^ for 
w — )■ oo. We could not find simple expression for C: it may depend on the initial data in a 
complicated way. 

Combining these results with the late time behaviour result for i^o 7^ in Eq. ( I29ll . the 
following expression for the late time behaviour of 0|r=M is consistent with all of our numerical 
results: 

2.00iJo C , , 

(p\r=M ~ h ^ , V OO . (30) 

All of our initial data has vanishing NP constant. We will explain below why the coefficient 
of f ~^ must be modified for initial data with non-vanishing NP constant. 

^ The exponents do not change up to the third decimal place even if we change the grid size as ((5C/, bv) — )■ 
(4(5C/, 5(Sw). Thus, the numerical error mainly comes from truncations. The truncation error can be estimated 
as OiXjv) ^ 10^"^. The exponents obtained here are reliable up to around second decimal place. 
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Figure 5: We depict (j)\r=M, dl(f)\r=M and d^(f)\r=M as functions of v for / = and Hq = 0. 
We take several initial data: a = 0.05,0.1,0.15 and /i = for an outgoing wavepacket and 
{a',fi') = (3.0, 10.0) for an ingoing wavepacket. For the ingoing wavepacket there is a period 
of damped oscillation (quasinormal ringing) before the power-law tail behaviour takes over. In 
all cases we find that, for large v, 4>\r=M decays as v~^, while d^(l)\r=M approaches a non-zero 
constant and d^(j)\r=M blows up as v. This implies that there is an instability even for initial 
data with -ffo[V'] = 0. 



14 




Figure 6: Time dependence of (p\r=M and df(j)\r=M for / = 1 and Hi ^ 0. We consider outgoing 
wavepackets with (a, /i) = (0.1, 0), (0.1, —0.05), (0.05, 0). After some damped oscillations 
(quasinormal ringing), 0|r=M exhibits power-law decay while 9^0 blows up linearly in v. 

3.4 Numerical simulation for / = 1 

3.4.1 Non-zero Aretakis constant 

Now we consider the case of the / = 1 partial wave, first with H\ ^ ^. For the outgoing 
wave we choose the parameters in the initial data to be (cr, /i) = (0.1, 0), (0.1, —0.05), (0.05, 0). 
In Fig. [6] we plot (^\r=M and 9^<^|r=A/ as functions of v. We can see that 0|r=M decays and 
9^(/>|r=A/ blows up linearly in v, as proved by Aretakis. 

Next we consider the late time behaviour of 0|r=M- Fig. [6] is consistent with a power law 
decay of </'|r=M of the form f " for large v. Fitting the absolute value of (^\t=m to such a power 
law in the range 1000 < v < 2000, we obtain the exponents a = —1.993, —2.001, —1.995 
for (cr, /i) = (0.1,0), (0.1,-0.05), (0.05,0), respectively. These results suggest that the scalar 
field decays as v""^ along the future horizon. To determine the coefficient of the power law 
decay, we fit the 4>\r=M to the function aHi/v"^ + l3/v^ for 1000 < v < 2000. Then we find 
a = 0.665, 0.665, 0.661 for {a, n) = (0.1,0), (0.1, -0.05), (0.05,0), respectively. Hence, for all 
cases that we have considered, the late time behaviour of (f)\r=M is given by 

(P\t=m 5 — V oo . (31) 

We will see in section S] that Eq. flM]) can be obtained from an AdS2 calculation. In Ref . |1] 
it was shown analytically that the / = 1 mode of the scalar field decays at least as fast as 
^-3/4 ^Qj, large v\ this is consistent with our numerical result. 

3.4.2 Zero Aretakis constant 

Now consider the case / = 1 with i^i = 0. For the outgoing wavepacket we choose the 
parameters a = 0.05,0.1 and /i = a{a — V(t2 + 4M2)/(2M) to set Hi = 0. For the mgomg 
wave we choose the parameters (a', /i') = (3.0, 10.0). We plot the time evolution of (p\r=M and 
9^0|r=M in Fig. [3 This is consistent with the following behaviour: 0|r=Af decays and (9^0|r-=M 
blows up linearly in v. Thus, we can conclude that there is an instability for / = 1 even when 
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Figure 7: Time evolution of 0|r=M and d^(j)\r=M for various initial data with / = 1 and 
Hi = 0: outgoing wavepackets with a = 0.05,0.1, fi = a{a — ya^^+^M^) / {2M) and an 
ingoing wavepacket with (o"',/i') = (3.0, 10.0). We find that 0|r=A/ decays and 9^0|r=A/ blows 
up linearly in v. 

the Aretakis constant vanishes. However, just as for / = 0, it affects quantities with one more 
r-derivative than in the case with non- vanishing Aretakis constant. 

Now consider the late time behaviour of (j)\r=M- Fitting the absolute value of <^|r-=M to a 
power law decay w", we obtain the exponents a = —3.001, —2.990 for the outgoing wave with 
cr = 0.05, 0.1 and a = —2.990 for the ingoing wave, respectively. This suggests (f)\r=M ~ C/v^ 
as f — oo for some constant C depending on the initial data. 

Combining this late time behaviour with the result for //i 7^ in ( 13T]) . the following late 
time behaviour for / = 1 initial data is consistent with all of our numerical results: 

0.66ifi C 

r=M 5 \- ^ V ^ 00 . (32) 

As we will explain below, we expect the coefficient of to be modified for initial data with 
non- vanishing NP constant. 

3.5 Late time tails outside the horizon 

Now we consider the late time behaviour of the scalar field outside the horizon. We 
consider outgoing wavepackets with (a, /i) = (0.1, —0.1), (0.1,0), (0.05, —0.05) for / = 0, 1,2, 
respectively. In Fig. |S]we plot the time dependence of at r = 1.5M. Initially the field is 
supported near the horizon (r = M) so the amplitude at r = 1.5M is small. The amplitude 
grows as the field disperses outwards from the black hole, and then decays as the field disperses 
to infinity. 

Our results are consistent with power law decay at late time. Fitting to f " we obtain the 
exponents a = —2.008, —4.006, —6.026 for / = 0,1,2, respectively. For the fitting, we use 
numerical data in the range 1000 < v < 2000 for / = 0, 1. For / = 2, we used data in the range 
800 < V < 1000 since the numerical calculation breaks down aX v > 1000. (We observed small 
oscillations which depend on the resolution.) 
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Figure 8: Time dependence of scalar field (f) outside the horizon: r = 1.5M. We consider 
outgoing wavepackets with {a, fi) = (0.1, —0.1), (0.1, 0), (0.05, —0.05) for / = 0,1,2, respec- 
tively. For such initial data, the conserved quantities are non-zero. After an initial period of 
damped oscillation (quasinormal ringing) we find power law decay 01^=1 5^/ ~ Ct'~^'~^. 

From these result, we extrapolate that the late time behaviour of the Ith partial wave of 
the scalar field outside the horizon is given by 



where rg > M for some constant C depending on tq and the initial data. We will reproduce 
this power-law from an AdS2 calculation in the next section]^ 

Ref. [16] argued that the decay outside the horizon should be proportional to t;~(2«+At+i) 
where = 1 if the NP constant is non-zero and /i = 2 otherwisejf] Our results show that this 
is not quite correct: instead we should have /i = 1 if either the Aretakis constant or the NP 
constant is non-zero (Ref. [16] considered only initial data with vanishing Aretakis constant.) 
This makes sense given the existence of the conformal isometry discussed in section 12. 3^ which 
interchanges a solution with non-vanishing Aretakis constant with one with non-vanishing NP 
constant. 

Ref. [16] found that the decay along the horizon is proportional to for vanishing 

Aretakis constant but non- vanishing NP constant. This implies that, in general, the coefficient 
of in (130 p and in (152]) will contain an extra term proportional to the NP constant 
(which is zero for our initial data). The decay along the horizon is proportional to v~^''^^\ as 
found in Ref. [16], but with /x defined as above. 

The decay at fixed r outside the horizon involves a different power of v than the decay on the horizon. 
This is consistent: C will diverge as r — ^ tq. 

^ Ref. [H] said that /x = 1 for an "initially static moment" and ji ~ 2 otherwise. But, according to Ref. 
|18) . the precise notion of an "initially static moment" is a non- vanishing NP constant. 




Cv 



-21-2 



V 00 



(33) 
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4 Scalar fields in AdS2 



In this section we will show that certain features of the late time behaviour on and outside 
the horizon of extreme RN can be deduced by looking purely at the near-horizon geometry. 



4.1 Coordinates 

The near-horizon limit of extreme RN is AdS2 x S'^ where each factor has radius M. 
Taking M = 1, the AdS2 metric in static coordinates is 

ds^ = -r^e + ^ . (34) 

Defining 

u = t + -, v = t-- (35) 

we obtain the metric in ingoing Eddington-Finkelstein coordinates 

ds"^ = -r^dv^ + 2dvdr . (36) 

These coordinates are regular at the future Poincare horizon r = 0. We can also use {u,v) as 
coordinates: 

4 

ds^ = —-. rT^dudv . (37) 

[u — vY 

Now define Kruskal-like coordinates 

U = tan"^ u, V = tan"^ v (38) 

to obtain 

ds' = - . ^ ^ dUdV . (39) 
sm ({/ — V) 

These are global coordinates on AdS2- The Penrose diagram for AdS2 is depicted in Fig. |9l 
The region covered by the static coordinates is — 7r/2 < V < U <7r/2. The left timelike 
infinity is U = V + ir, the right timelike infinity is U = V. The future Poincare horizon of 
the static patch is U = 7r/2, the past Poincare horizon is = —tx/2. 



4.2 Massless scalar in AdS2 

Consider first a massless scalar ip field in AdS2- Since its equation of motion is conformally 
invariant, we can write down the general solution: 

i^ = f{U)+g{V) (40) 

for arbitrary functions /, g. To compare with our earlier results, we need to calculate the 
Aretakis conserved quantity. Writing out the wave equation in {v, r) coordinates and evalu- 
ating at r = (or taking the near- horizon limit of (^) shows that the following quantity is 
conserved along the future Poincare horizon 

Ho[^] = idr^)r=0 . (41) 
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Figure 9: Penrose diagram for AdS2. 

Let's evaluate this in {U, V) coordinates. Converting from {v, r) coordinates to {U, V) coor- 
dinates gives 

d_ ^ sm'jU-V) d 

dr 2 cos2 V dU ' ^ ' 

Hence 

^0 = -^/'(vr/2) . (43) 
Now consider the second r-derivative of the field at the horizon: 

{d^4^)r=o = ^r(vr/2)tany + ^r(7r/2) = -H^v + ir(7r/2) (44) 

so when i^o 7^ we have linear growth in v, just as for extreme RN. On the other hand if 
Hq = then we see that dlip\r=o is a constant; by taking another r-derivative we see that if 
this constant is non-zero then d!^ip\r=Q will blow up linearly in v, just as we found numerically 
for extreme RN. 

Is this an instability of the scalar field in AdS2'^ No: working in (f/, V) coordinates one 
sees that tp and all of its derivatives are bounded. The blow-up in [v, r) coordinates is simply 
a coordinate effect, it has no invariant significance. By contrast, in extreme RN, the {v,r) 
coordinates are a preferred set of coordinates that follow from the symmetries of the spacetime. 
The function r is defined using spherical symmetry and there is a unique Killing vector field 

that is timelike and canonically normalized {K'^ = —1) at infinity which is used to define 
V. In AdS2 there exist infinitely many timelike Killing vector fields analogous to K'^. Hence, 
unless one has some reason to regard one of these as preferred, then there is no invariant way 
of introducing {v,r) coordinates and hence no instability. 

Now consider the late time behaviour of the field. The value of the field on the Poincare 
horizon is 

^l;i7^/2,V) = fi7^/2) + giV). (45) 
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Converting to {v,r) coordinates, the late time behaviour of the field at the horizon {v — )■ oo, 
r = 0) is 

V;(vr/2, V) = /(7r/2) + ^(7r/2) - + ^ ( ^) ■ (46) 

So far we have imposed no boundary condition on %p. Note that the value of ip on the right 
timelike infinity iqj 

^^{U)^f{U)+g{U) (47) 

where we use U as the coordinate at right timelike infinity. We now impose the boundary 
condition that ipooiU) and 'ip'^iU) decay at late time at right timelike infinity: 

V'oo(f/)^0, VL(f^)^0, as U^7i/2, (48) 

which is equivalent to 

g{rr/2) = -/(7r/2), 9'{n/2) = -f{n/2) (49) 

and hence 

i;{n/2,V) = -^ + o(\] . (50) 

This agrees to good accuracy with our numerical results for extreme RN (summarized in 
equation fl5U|l )]1 The late time behaviour outside the horizon (w — )■ oo at fixed r > 0) is, if we 
make the extra assumption ip'^{U) as U 7i/2, 

= + O f^) ■ (51) 

This is in good agreement with ( 133|) . Note that our assumptions on ipoo are all satisfied if 
there exists Uq < ti/2 such that ifjoo vanishes for U > Uq. 



4.3 Massive scalar in 74^5*2 

Now consider a scalar field of mass m in 74^52. Since the general solution of the wave 
equation is not so simple in the massive case, we will proceed more heuristically than we did 
above. We assume that, in a neighbourhood of f/ = V = n /2 (i.e. late time at right timelike 
infinity) the field satisfies "normalizable" boundary conditions 

iIj = {U -V)^"^F{y) + 0{{U -V)^"^+^) (52) 

where FiV) is an arbitrary function and 



^n^^- + ^jm^ + -. (53) 
^ In the language of AdS/CFT, il^ooiU) is the source for the operator dual to ij). 

^Recall that this result applies only for the case of vanishing NP constant. If the NP constant is non-zero 
then there is radiation present near infinity in extreme RN. Maybe one could reproduce results for this case 
by modifying the boundary condition for ^/^^ in (|48| . 
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Taking the near-horizon hmit of the Ith multipole ipi of a massless scalar in extreme RN 
(setting M = 1) results in a scalar field of mass m in AdS2 where 



= /(/ + !), / = 0,1,2,... (54) 

and therefore = 1 + 1. In this case, the conserved quantity at the Poincare horizon is 
simply 

Hi[^lj] = (9t+V)r=o • (55) 

We can now calculate Hi[ip\ at large v: r = 0,t'->oois equivalent to U = 7i/2, V — t- 7r/2, 
and we can use the behaviour (l52l) of ip near the boundary to obtain 

Hence, near U = V = n/2 we hav^ 

^ = J^l^H, {-2{U - V)y-'' + 0{{U - Vy+') . (57) 

From this expression we can calculate the late time behaviour along the Poincare horizon in 
(w,r) coordinates (r = 0, w — )■ oo): 



For / = 1 this agrees to good accuracy with our numerical result (summarised in ( 132|) ) for the 
full extreme RN geometry. The late time behaviour outside the horizon (r > fixed, v — )■ oo): 

i! ,. / 4 ^ 



in agreement with f l33p . 

Why are these AdS2 results in such good agreement with our numerical results for the full 
extreme RN geometry? One might have thought that the late time behaviour at the horizon 
would exhibit some dependence on the geometry outside the horizon, which is different in 
74^5*2 and extreme RN. Maybe the numerical coefficients in fl5U]l and fl5^ would deviate from 
the AdS2 values if they were computed to higher accuracy. Or maybe all that matters at late 
time is the geometry exactly at the horizon, which would explain the agreement. It would be 
nice to understand this better. 

Observe that taking the near-horizon limit of a spherically symmetric scalar of mass m in 
extreme RN (again with M = 1) also results in a scalar of mass m in AdS2- From the above 
we deduce that if = n{n + 1) where n is a positive integer the above calculations remain 
valid (with the replacement I ^ n). In particular, we learn there are conserved quantities 



^ It follows that, in AdS/CFT, the conserved quantity Hi is proportional to the late time limit of the 
expectation value of the operator dual to '0. 
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for certain massive scalar fields in the near-horizon geometry of extreme RN. In Sec. 15.11 we 
show that in fact a massive scalar field in the full extreme RN also possesses such conserved 
quantities, for these specific values of the mass, and show they lead to instabilities. 

For general m it does not appear possible to construct conserved quantities on the horizon; 
however from (152|) we may still deduce the late time behaviour. As f — )■ oo we find that on 
the horizon 

for all k > 0: we see this blows up for k > [AmJ + 1, where [ J denotes the integer part. In 
Sec. [5] we will demonstrate numerically that instabilities of this nature indeed occur in the 
full extreme RN for certain m such that is not integer (despite the absence of a conserved 
quantity). For completeness, the late time behaviour outside the horizon at fixed r > is 

^ - -IZ^ ■ (61) 



5 Massive scalar field 

In this section we consider the evolution of a spherically symmetric massive scalar field in 
the extreme RN background. First we will show that, for special values of the mass of the 
scalar field, a tower of conserved quantities exists just as for the massless scalar field. These 
can be used to deduce the existence of instabilities. Then we will solve the equation of motion 
numerically as in the massless case to investigate the evolution in more detail. 



5.1 Conserved quantities on the horizon 

We start by writing down the massive wave equation f l2^ for a spherically symmetric 
massive scalar field in ingoing EF coordinates 

2r9„a,(r^) + dr (Aa,^) - mV^^ = . (62) 

We find that the construction of the conserved quantities for the massless scalar field (fTOj) can 
be generalised to a massive scalar for special values of the mass: 

= n{n + 1)M-^ (63) 

where n is a positive integer. This gives a very similar term in the equation of motion to 
the centrifugal barrier of a massless scalar field, so the construction of conserved quantities 
proceeds analogously to that of the massless scalar with n ■v^ I. Explicitly, we find 

H4i^] ^ i- {d: [fn{r)dr{rmr=M (64) 
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is conserved along H"*", where is a polynomial defined bj|§ 



/c=0 



(65) 



whose coefficients are recursively determined by cq = l,ci 



n and for 2 < k < n 



n{n + 1) [2{n -k + l)cfc_i + (n 



/i;)Cfc-2] 



(66) 



k{n-k + 2){2n-k + l) 



It follows that generically not all of the quantities d^ip for < j < n + 1 can decay at late time 
on 'H'^. By analogy with the / > massless scalar, we expect, and our numerics conffims, 
that d^lip decays at late time on 7/"*" for j < n. Then d^^^ip does not decay and d^^'^'^ip blows 
up linearly along "H^. 

It is interesting to ask what happens for general mass m. Although we do not find any 
conservation laws in this case, by continuity in the mass parameter one might expect that, 
generically, d^^'^^ip does not decay, and d'^^'^ip blows up, at late time on "H^ where n is given 



where is given by Eq. f l53|) with m replaced by Mm and [ J represents the integer part. 
We confirm this expectation numerically in the next section. 

5.2 Numerical results 

Consider the massive scalar wave equation in the coordinates developed earlier fl25l) . Using 
the numerical method introduced for the massless scalar field discussed in section [3|, we study 
the time dependence of a massive scalar field. 

5.2.1 Non-zero conserved quantity 

In this section, we study the time evolution for the case corresponding to that studied in 
section 13.3. for which the conserved quantities 0641) are non-zero. 

First of all, we show the time evolution of (j){v,r) = ripiy^r) on and outside the horizon 
for n = 1 {{MmY = 2) using the outgoing wave initial data fl27|l with (ex, /i) = (0.1, —0.1). In 
Fig. [TOl we plot 0(f , r) and the quantity Hi{v, r) defined by 



which coincides with the conserved quantity ( 164|) for n = 1 at the horizon. Panel (a) shows 
that 0(f,r) quickly decays on and outside the horizon. Panel (b) shows that Hi{v,r) stays 
constant on the horizon while it decays to zero outside the horizon. This property is consistent 
with d^(j)\r=M being constant and df.(f)\r=M diverging linearly in v at late time. 

^In fact any smooth /„(r) whose derivatives /n'^^|r=A/! for /c = 0, . . . ,n, coincide with those of this poly- 
nomial will do. 



by 



n = 



[Am\ - 1 



(67) 




(68) 
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Figure 10: (j){v,r) and Hi{v,r) for the massive scalar field with {MmY = 2. 



In Fig.[TTl we show 4'\r=M and d'^(f)\r=M for {MmY = 2 for various initial data with Hi ^ 0. 
To generate this plot, we used the outgoing initial data f E7|) . Panel (b) shows that dl(t)\r=M is 
linearly divergent with respect to time, as expected from existence of the conserved quantity 
Hi (this could be proved analytically in the same way as for the massless scalar with / = 1). 

In Panel (a), we find that, at late time, shows oscillation with frequency u = m. 

Oscillations with this frequency have been seen in previous studies of massive scalar fields in 
black hole backgrounds [2D] . It is convenient to separate the value of the field at the horizon 
into its mean value and an oscillating component: 

(P\r=M{v) = <I>~{V) + r^'-ivt^M ' (69) 

where (f)\^l^{v) = J^^^j^ (j)\r=M{v')dv' with T = 2'ji/m. At late time, the oscillating 
component decays more slowly than the mean value. On the other hand, we observed that 
for d^-^(f)\r=M, the oscillations are subdominant at late time for generic initial data. This 
suggests that the instability is associated to 0™'''^°. 

Next we examine in more detail the behaviour of each component in (j)\r=M- First, we 
show the behaviour of ) for {Mmy = 2 in Fig. [121 We find that ) decays by 

a power law, and its value is governed by the conserved quantity Hi as 



I mean 

r\-/ 



OMHi 



M T" 



(70) 



for several different choices of outgoing wave initial data. More precisely, we calculated 
the late-time behaviour of 0|r=A/ for outgoing wave initial data for (ex, /x) = (0.1,-0.1), 
(0.15,-0.1), (0.05,-0.1), (0.1,-0.15) and (0.1,-0.05). We fitted to the function 

aHi/v"^ + for 600 < < 1000, and found a = 0.664, 0.665, 0.663, 0.663 and 0.663 
respectively for these initial data. Note that the behaviour of Eq. (170]) was seen also in the 
case of massless scalar field with I = 1 (See Eq. (13T]) ). 

Next, we examine the behaviour of (f)\°^^\,j. Since this component oscillates with frequency 
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Figure 11: Time evolution of 4'\r=M and d^(j)\r=M at the horizon for (Mm) 
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Figure 12: \v (j)]^^'^ / Hi\ for outgoing wave initial data for (cr, /i) 
(0.05,-0.1), (0.1,-0.15) and (0.1,-0.05). For any initial data, \v 
at late time. 
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m at late time, it is useful to focus on its amplitude y4|°!f]y-(w) defined by 



A\ 



r=M 



(71) 



In Panel (a) of Fig. [131 "we plot A\"!^"i^j{v) in the case of outgoing wave initial data with 
= (0.3,-0.1) for {Mmf = 2,4 and with (0.15,-0.1) for {Mmf = 2. Fitting the 
data to the function Cf", we find the power of decay is a = —0.84, 0.83 and 0.82 for our 
different initial data. These results are consistent with the scaling A\°^";^^{v) oc v~^^^. We plot 
also the oscillation amplitude of (j){v,r) outside the horizon in Panel (b) for the same initial 
data, which shows the same scaling A\°^"^ smI"^) ^ More precisely, fitting the late time 

behaviour with Cv"" gives a = —0.84, 0.84 and 0.82 for our different initial data. Note that 
y~5/6 gQaiing for the late-time tail of a massive scalar field outside a black hole horizon has 
been discussed previously in Refs. [201 EI] • The constant C in these fits does not appear to 
be simply related to the conserved quantity Hi. 
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Figure 13: Panel (a): Amplitude of ^^^1 A\°^"f^{v), in the case of outgoing wave initial 

data with (a, /i) = (0.3, -0.1) for {Mmf = 2,4 and that with (0.15, -0.11) for (Mm)^ = 2. 
Panel (b): A\°^"{v) at r = 1.5M for the same initial data. Numerical results suggest that 
A\°^"]^(v) and A\°^f[ 5a/(^) become proportional to v~^^^ at late time for any initial data. 

We also comment on the dependence of 4>\r=M on the scalar field mass m. For this purpose, 
it is useful to focus on time dependence of Hq and Hi, which are conserved when [MmY = 
and 2, respectively. We summarize their time dependence in Fig. [TH These quantities are 
dominated by 0™*^'^°^ and show power law decay or growth at late time. Assuming the power- 
law dependence Hn oc v^", we can obtain the powers Pn by fitting the numerical data. We 
summarize the resulting p„ in Fig. [151 The numerical results strongly suggest for general m 

Hn^v^% p„ = n + 1 - A„ . (72) 

Using the results in section 14.31 we can see that this is the same scaling as appears for a 
massive scalar field in an AdS2 background. 
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Note in Panel (b) of Fig. [TUthat pi is positive for [MmY = 1 which confirms the existence 
of an instabihty in this case, for which the analytic argument above does not apply. The 
instability involves blow-up 9^0|r=M5 in agreement with the prediction of f E7|) . 
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Figure 14: dependence of Hq{v) and Hi{v) for [Mm) 



0,2, 




Figure 15: dependence of p„, where we assume if„ oc w^". In both panels, the 

numerically obtained p„ for (Mm)^ = 0,1,..., 8 are shown with analytic curves Pn = 
ri + (l — a/1 + 4(Mm)^) /2. For any data point, relative error of the numerical value from 
the analytic value is ~ 5 x 10~^ at most. 

We close this section with some comments on the behaviour of derivatives of the field 
outside the horizon. In Fig. [161 we show d^(f) on the horizon for (Mm)'^ = 2, and also that 
quantity off the horizon for {MmY = 2, 8. As expected from the general argument, 
converges into a constant at late time. On the other hand, 9^0 outside the horizon shows 
damped oscillation whose amplitude decays as v~^^^. This behaviour appears universal for 
any scalar field mass m and also order of the derivative, that is, the envelope of the damped 
oscillation of 9^0 is proportional to v~^^^ for any k and m. Conservation or blow up of 9^0 
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is a property that occurs only on the horizon of an extreme black hole, and away from the 
horizon these quantities simply decay (as proved by Aretakis in the massless case). 




Figure 16: \d^(j)\ on the horizon for {MmY' = 2 contrasted with that at r = 1.5M for 
{MmY = 2,8. for (Mmy = 2 converges to a constant at late time as suggested by 

the general argument, while 9^0 off the horizon decays as oc t>~^/^ irrespective of the mass m. 



5.2.2 Zero conserved quantity 

In this section, we briefly comment on the cases Hn[ip] = for the massive scalar field. 
We focus on the case n = 1 ((Mm)^ = 2), but the results for general n are analogous. 

First, we show the time dependence of (l>\r=M in Fig. [T71 To set ifi = 0, we used the 
outgoing wave initial data with = —a [a + Va^TUP) /(2M). The behaviour for the 
outgoing wave initial data is similar to the Hn 7^ case studied above in that 0™^'^'^ and (f)°^'^^ 
are non-negligible at r = M. On the other hand, (f)\r=M{v) for the ingoing wave initial data 
is dominated by the oscillatory component with negligible contribution from 0™°^'^. 

We summarize the behaviour of the mean and oscillatory components of 4>\r=M in Fig. [181 
First, we find from Panel (a) that the mean component shows the power-law decay 01^?^^ ^ 
v~^, which is one power faster than the Hi ^ case studied above. This property is similar to 
that of the massless case with / = 1, and we can show that d^(f)\r=M diverges linearly in time 
(as it did in the massless case). Hence, just as in the massless case, vanishing Hi requires one 
more r-derivative to see the instability. 

Next, we find from Panel (b) that the amplitude of the oscillatory component shows the 
power-law decay A\'^^^\,j{v) oc at late time. This is the same power law as the Hi 

case shown in Fig. [131 and thus the time dependence of Aosd seems to be insensitive to the 
value of the conserved quantity. 

To sum up, the behaviour of the component 0|^^ of a massive scalar field with (Mm)^ = 2 
can be described by Eq. ( l32l) . which was proposed for the massless scalar field with / = 1. For 
general value of the scalar field mass m, our results for 0|^^/ are consistent with a power-law 
decay v~^"^ for generic initial data, or y-^^^i-i ^qj. initial data with Hi = 0, where is 
defined below Eq. (!72|) . On the other hand, the amplitude of the damped oscillation of 4>\r=M 
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100 

Figure 17: (f)\r^]x.j{v) for [Mrnf' = 2 with Hi = 0. We used the outgoing wave initial data 
with a = 0.05 and 0.1, and also the ingoing wave initial data with (a, fi) = (3.0, 10.0). For the 
outgoing wave initial data, n is fixed as /x = —a (a + \/a'^ + / (2M) by the requirement 

of Hi = 0. 

or the derivative d^(j) outside the horizon, appears to be proportional, at late time, to v~^^^ 
irrespective of m and also of order of the derivative. 

6 Gravitational and electromagnetic perturbations 

In this section we show that an analogue of Aretakis' scalar field instability occurs for 
the coupled linearised gravitational and electromagnetic perturbations of the extreme RN 
solution. To demonstrate this we will employ past work on the linearized perturbations of 
extreme RN. Several formalisms were developed for studying this problem, see [221 1231 [211 
[25l [26l [271 [28] and also |29] which relates the various methods. These works established the 
mode stability of the RN black hole (including the extreme case). 

We will use Moncrief 's gauge invariant formalism [23l [2H [25] due to the remarkable sim- 
plicity of the resulting perturbation equations. This approach is based on the Hamiltonian 
formulation of the Einstein-Maxwell equations, where the perturbation variables are certain 
gauge invariant combinations (under both infinitesimal diffeomorphisms and electromagnetic 
gauge transformations) of components of the metric and gauge field. We will first briefiy ex- 
plain Moncrief 's formalism and variables, then we present conservation laws for these variables 
and associated instability results. 

6.1 Moncrief 's gauge invariant formalism 

Moncrief simplified the perturbation equations for the RN solution by employing a reduc- 
tion method based on the Hamiltonian formulation of Einstein-Maxwell theory of ADM. 

First recall that in the Hamiltonian formulation one chooses a preferred time function t 
and decomposes the spacetime fields on hypersurfaces of constant t, with coordinates 
for i = 1,2, 3, as follows. The metric g^i, is written in terms of the lapse = —(5'**)^''^, shift 




29 




1e-05 ^ CT=0.05 

CT=0.10 

1e-06 1^ ■ — 

10 



N 



0.0001 



50 



100 



500 



10 



100 



1000 



V 



V 



(a) 01 



mean 



(b) A\ 



OSCl 



r=M 



Figure 18: Panel (a): ) of {MmY = 2 and i/i = for the outgoing wave initial data 

with (J = 0.05 and 0.1. Fitting the curves with the function Cv"", we find the power of the 
decay to be a = —3.01 and —2.99 for each initial data. We do not show the result for the 
ingoing wave initial data since (pl'^^M is almost vanishing compared to 4'\r=M it. Panel (b): 
Amplitude of the damped oscillation of (f)\°!^^\.j{v), A\°^^\,j{v), for the ingoing wave initial data 
with (cr, /i) = (3.0, 10.0) and outgoing wave initial data with a = 0.30,0.15. Fitting the data, 
we find the decay —0.85, —0.84 and —0.84. These results suggests that A\°^^\.j{v) becomes 
proportional to v~^^^ at late time, which is the same power law as in the Hi case studied 
above. 

Ni = Qti and the induced 3-metric Qij on S^, whereas the gauge field is decomposed in 
terms of the scalar potential Aq and the spatial vector potential Ai. The canonically conjugate 
momenta to Qij and Ai are the tensor densities vr*-' = ^{K^^ — Kg^^) and = N ^F^'^ 
respectively, where Kij is the extrinsic curvature on S^. Variation of the action functional with 
respect to these canonically conjugate pairs leads to dynamical evolution equations, whereas 
variation with respect to A^, Ni, Aq leads to the constraint equations Ti = Hi = £ = 0. 

We are now ready to give a brief review Moncrief's reduction method, see [25] for more 
details. The perturbation equations can be derived by taking the second variation of the 
action functional evaluated on an exact solution. The canonical variables for the linearised 
ADM equations are simply the first order perturbations of the canonical variables {gij,7r'^^) 
and {Ai,E^), which we will denote by {6gij,6'K^^) and {6Ai,6E^). Similarly, the constraint 
equations for the perturbed problem are the linearised constraints 671 = 6W = 68 = of the 
exact problem; these must be conserved in time as a consequence of the evolution equations for 
the canonical variables. There are two crucial facts about the linearised constraint functions 
61-1, 6'W, 6£: they generate infinitesimal diffeomorphisms and electromagnetic gauge transfor- 
mations, and mutually commute under the Poisson bracket. These properties allow one to 
perform a canonical transformation to new canonical variables such that 6%, 6W, 68 are a sub- 
set of the momenta and their conjugate variables are cyclic and gauge- dependent. Hence the 
remaining canonical pairs, which by construction now will commute with 671, 6'W', 68, must 
be gauge-invariant. In terms of these new canonical variables the problem greatly simpli- 
fies: the evolution equations for the gauge-invariant pairs decouple from all gauge dependent 
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quantities and the constrains can be eliminated. 

Moncrief applied the above general method to perturbations of the RN solution: naturally 
the time function t was chosen to correspond to the static time coordinate. The background 
data is thus (we don't yet assume extremality) 

= F{r) = l- — + ^ , g,,dx'dx^ = ^ + r^dn\ (73) 

iVi = 7r^^ = 0, Ai = 0, E^ = 2gsin^. (74) 

Spherical symmetry of the background solution allows one to further simplify the problem. 
One expands all perturbation variables in terms of spherical harmonics Yjm which then leads 
to decoupled equations for each value of / = 0, 1, 2, . . . . We now describe the basis used to do 
this, which is essentially that used by Regge and Wheeler. 

There are three classes of perturbation depending on whether they transform as a scalar, 
vector or tensor on S*^. A basis for each of these (written covariantly) is given by the following 
harmonics: 

Scalar (/ > 0) : Y = Yim where - = /(/ + l)Y (75) 
Vector (/>1): = -^^.F, fa = daY (76) 
Tensor (/ > 2) : e.^ = e^,'Vb)V,Y, U = V.VfeF, = Y{dn^)ab, (77) 

where dfl^, Vq, iab are the round metric, the associated connection and volume form on the unit 
S"^ and a,b, . . . are S"^ indices. As always the Yim are orthonormal with respect to the unit 
measure; from the properties of Y one may deduce corresponding properties for the vector and 
tensor harmonics. These harmonics all have a well defined parity under the discrete isometry 
which reverses the orientation. In particular Y, fa, fab, dab have parity (— 1)^ whereas ea,eab 
have parity (— l)'"*"^. Traditionally these are referred to as even or odd parity respectively; 
for each / the perturbation equations for the two types of parity therefore decouple. We note 
that since we are in four dimensions perturbations of the tensor components of the Maxwell 
two form can be expanded in the basis Yiab, which is an odd parity harmonic. 

By Birkoff 's theorem any / = perturbation must correspond to the RN family of solutions. 
Perturbations which correspond to a slowly rotating Kerr-Newman black hole are contained 
in the / = 1 odd type, as discussed below. Hence any instability must either be an / > 2 
perturbation or a non-stationary I = 1 perturbation. 

Moncrief showed [231 1211 ISSj that each type of perturbation can be fully described by 
either one gauge invariant variable (/ = 1) or a pair of gauge invariant variables (/ > 1). 
We describe these below and more fully in Appendix [Bl Remarkably, each of these variables 
satisfies a decoupled wave equation of the general form: 

2 

- ^a^V + dr{r-'^Adrij) = r-^W{r)ij (78) 

where A = r^F, for some explicitly known function W{r). However, one shortcoming of 
Moncrief's approach is that since the various are defined using static coordinates, it is 
not obvious if they are smooth on the horizon (since those coordinates breakdown there). 
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In Appendix [B] we show that in fact all the relevant gauge invariant variables if) are indeed 
smooth on the future horizon H"*". Hence we may convert to EF coordinates resulting in the 
master equation 

2d^dr^ + drir-^/^dr^) = r-^W{r)'iP . (79) 

All terms in this equation are smooth on "H"^ which is at the largest root r = r+ of A. 

The odd parity perturbations are simplest, so let us describe the gauge invariant variables 
in this case (see Appendix [B] for more detail on their construction). 

General / > 1 odd perturbations of the gauge field must be of the form 6Aa = Aca for some 
function A. It turns out ttj = A is one of Moncrief's gauge-invariant canonical momentum 
variables. This variable can be defined invariantly as follows: 

where S'^{v,r) a sphere of constant {v,r) in EF coordinates and the second equality follows 
from integration by parts. Under a gauge transformation 6F = L^F = di^F where is a 
vector generating infinitesimal diffeomorphisms. For odd perturbations C./idx'^ = C{v,r)eadx°' 
and it is readily checked that this implies 6 Fab = 0. Hence vr/ is indeed gauge invariant and 
the above definition shows it is also smooth on the horizon. 

The relevant variable to describe / > 1 odd perturbations of the metric turns out to be 
Pi, which is the conjugate momentum to the metric perturbation 6gra = hiCa- In Appendix 
[B] we give details of general odd perturbations of the metric and show that 

Pi= I dn in* (djgra - dr6g,a + , (81) 

where here 5gfj,iy is written in EF coordinates, which explicitly shows pi is a smooth quantity on 
Ti^. It is easily checked this is gauge-invariant directly. For odd perturbations an infinitesimal 
diffeomorphism must be of the form Sg^u = V{fj.^u) where again C,^dx^ = C{v,r)eadx°'. Hence 
Sgva = dyC and 6gra = drC — 2C/r; it follows that 5tii = as claimed. 

Moncrief showed that / > 1 odd perturbations reduce to two decoupled wave equations 
fl79p for ip = P± where P± are certain constant linear combinations of ttj and 

TTg = rpi TTf . (82) 

From above we deduce that TTfjUg and hence P± are indeed smooth functions on the horizon. 

For / = 1 odd perturbations Moncrief showed that the quantity 6a = —-^{^r^Pi — '^Q'^f), 
which from the above is gauge invariant and smooth on H^, is in fact a constant. In this case 
the perturbation equations reduce to a single wave equation of the general form fl79|) with 



^ = Pf^^f~ w = 2{l + ^] . (83) 



From the above we deduce that Pf is a gauge invariant quantity which is smooth on the 
horizon (since both Tif and 6a are). Note that Pf = corresponds to a perturbation to a 
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slowly rotating Kerr-Newman black hole with rotation parameter 5a (explaining the notation); 
indeed this is the most general regular stationary solution in this case [2^ . 

The even parity perturbations are more complicated to describe. In Appendix |B] we 
explain their definition and show that the gauge invariant quantities appearing in Moncrief 's 
final decoupled equations (179|) are also all smooth on the horizon. We establish this by showing 
how to express them as linear combinations of the metric and field strength perturbations in 
EF coordinates. 



6.2 Conservation laws and blow up on the horizon 

As explained in the previous section, Moncrief 's perturbation equations can all be reduced 
to the form f l7^ for some gauge invariant function which is regular on T-C^ . We will now 
restrict to the extreme RN solution Q = M so A = (r — M)^ and r+ = M. Using this 
general form we now derive necessary and sufficient conditions for the existence of certain 
conserved quantities on the horizon. As we will show below, remarkably, these conditions are 
satisfied for all types oil >1 perturbation. Hence coupled gravitational and electromagnetic 
perturbations of an extreme RN solution possess a tower of conservation laws on the horizon. 
In turn, these conserved quantities lead to blow up results and hence instabilities for all types 
of perturbation. 

Let /(r) be a smooth function which is non-vanishing on the horizon and without loss of 
generality set f\r=M = 1- Multiply f l7^ by r'^f{r) and differentiate p > 2 times with respect 
to r and then set r = M. This gives an equation of the form 

p 

2d,[dP{r^fdrmr=M + Yl «fc^r'^lr-=M = , (84) 

fc=0 

where in particular 

ao=p{p + l)-W\r=M. (85) 

Hence the quantity in square brackets in fl8^ is conserved on the horizon if and only if = 
for all < k < p. In turn, this is equivalent to the following set of equations: 



/ 



W\r=M = P{P+1) (86) 
1 



(fc) 



2p + l-k 



L q=l 



7) 



where 1 < k < p. For each k > 1 equation (!87|) determines f^''^\r=M in terms of the lower 
order derivatives f^''~^^r=M,---,f\r=M- It follows that for all 1 < k < p the f^''^\r=M are 
determined (recall we have fixed f\r=M = !)• Hence by choosing / to be any function whose 
derivatives satisfy these constraints we deduce that 

Hpm^^mr'fdM\r=M (88) 

is conserved along if and only if (l86l) is satisfied. Notice that Hp depends only on the first 
p derivatives of / at r = M and is therefore independent of any specific choice of /. 
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This conservation law implies generic non-decay at late time on of a certain linear 
combination of ij) and its first (p+ 1) r-derivatives. Since we expect all derivatives of the field 
to decay outside^ "H"^, it follows that certain quantities with one more r-derivative will blow 
up at late time on "H"*". 

To see this explicitly, suppose we have a conserved quantity Hp[^] for some p > 2 as above. 
Now consider ( 18^ with p — p + 1 so 



P+i 



r=M 



^9) 



k=l 



where we have used ( 18 6 p to simplify the coefficient of df.~^^ip\r=M which in this case is non 
vanishing. Now (motivated by the results for the scalar field) suppose that ^^^|:\r=M — as 



f — )■ oo for < /c < p. Then flHHj) implies df^^ipl 



r=M 



Hence f l89|) implies 



r=M 



(90) 



as w — oo. Therefore if the constant Hp ^ 0, as will be the case for generic initial data, we 
deduce the linear blow up 



as u — )■ OO. Iterating this argument shows 



M2 



(92) 



for some non-zero constants Cp^n and all n > 1. 



Let us now apply these results to Moncrief 's perturbation equations. As we show below, 
remarkably, for all types of perturbation there exists some integer p > 2 such that ( l86l) is 
satisfied, as summarised in table [H 





/ = 1 odd 


/ > 1 odd 


I = 1 even 


/ > 1 even 




Pf 


P± 


H 


R± 


W\r=M 


6 


/(/ + !) + ! ±(2/ + 1) 


6 


/(/ + !) + 1± (2/ + 1) 


P 


2 


l±l 


2 


l±l 



Table 1: Conserved quantities Hp[%lj\ for Moncrief 's perturbations. 



Hence for all types of perturbation there exists a conserved quantity flHHj) and if this quantity is 
non-zero, there exists an associated blow up result for the corresponding transverse derivatives. 
Note that the smallest value of p, i.e., the fewest transverse derivatives, is p = 1 which occurs 
for / = 2. 

It would be surprising (and very interesting) if this expectation were false because then metric-Maxwell 
perturbations would exhibit behaviour that is worse than that of scalar field perturbations. The case / = 1 
is an exception for which the exterior solution might settle down to a perturbation within the Kerr-Newman 
family. 
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Consider the / = 1 odd perturbation equation fl83l) for which ip = Pf and iy|r=A/ = 6. 
In this case there exists a conserved quantity for p = 2: exphcitly one can use fl57|) to get 
f'\r=M = 3/2 and f"\r=M = and deduce that 

H2[Pf] = (d'^Pf + l^dlPf + ^,drP^ ^ . (93) 

For a generic initial perturbation, i?2[-P/] 7^ and hence drPf, d^Pf and d'^Pf cannot all decay 
at late time on T-L'^. In view of the results for a scalar field, the most likely behaviour is that 
Pf, drPf, d^Pf all decay along the horizon and hence d^Pf does not decay. In this case, from 
the above argument, d^Pj will blow up linearly in v on T-L^ . Hence we must have an instability. 
This is readily translated to a statement about explicit components of the Maxwell field: if we 
restrict to perturbations for which (5a = (e.g. by demanding vanishing angular momentum 
perturbation), then we deduce the best possible generic behaviour is that the 5*^ components 
5Fah, drSFab, d'^6Fab all decay, d^SFab does not decay and df6Fab must blow up linearly in v 
along "H"*". 

Now consider / > 1 odd perturbations. In this case the variables ip = P± are related to 
the gauge- invariant variables vrj,??^ defined in equations (ISUi) . (15^ and (15T1) . by 



1 



^2(2/ + !)(/-!)(/ + 2) 
/(/ + !) 



V2/ + 4 P: - V2/-2 P_ 



V2I - 2 P+ + V2/ + 4 P_ 



v/2(2/ + 1) 

whereas the associated functions appearing in the wave equation fl79|) are 



(94) 
(95) 



, 3M 4M2 (2/ + l)M , , 

W± = l{l + 1) + ± ^ ^ . 96 

Therefore W±\r=M = /(/ + 1) + 1±(2/ + 1) and hence W+\r=M = (/ + l)(/ + 2) and W. = /(/-I). 
We deduce that for P+ there exists a conserved quantity for p = / + 1 and for P_ a conserved 
quantity for p = / — 

For definiteness consider the / = 2 odd case so then we have conserved quantities: 

H^[P_] = (d'^P. + ^^drP.) (97) 

/ r=M 



M 

H,[P^] = (9,^P+ + ...),=M (9^ 



The recurrence relation (|87| can be solved explicitly in this case. For 

(,) 3fc!(-l)fc(/-fc + l)(2^-fc + 2) 

'"-^^ Mfci(/ + l)(2Z + l) 

and for P_ 

(,.) _ (-l)'=fc!(fc + l)(/-fc~l) 
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where the elhpsis denotes lower order r-derivatives of P_|_, whose exphcit coefficients are unim- 
portnt. 

For a generic initial perturbation these constants will be non-zero (and independent). 
Hence drP- and d^P- cannot both decay at late time on "H"*". One can relate P_ directly to 
the metric and Maxwell field perturbation by writing it as a linear combination of vrj and Hg 
using ( 194|) and ( p5|) . We can then use (182|) to write P_ in terms of vtj and pi. Recall pi is 
the gauge invariant combination of metric perturbation components (ISTj) and vtj is the gauge 
invariant Maxwell field perturbation f l80p . Hence we deduce that a certain gauge-invariant 
linear combination of the first two r-derivatives of these quantities generically does not decay 
at late time on "H+j^ 

As discussed above, it seems very likely that, for / = 2, the gauge invariant quantities and 
all of their derivatives will decay at late time outside H^. Hence, since we have non-decay of 
a certain quantity on H^, there will be a quantity with one more r-derivative which blows 
up at late time on "H^. The most likely scenario is that P_ and drP- will decay at late time 
along the horizon, in which case d^P- cannot decay and d'^P- will blow up at late time on 
"H^. Similarly, if P+, drP+, 9^P+, dfP+ all decay along "H^, then (9^P+ cannot decay and d^P+ 
will blow up at late time on "H^. It then follows that vtj, -Tr^, drTif, d^T^g decay, whereas 

a,^.,|„__^ (99) 
in which case from equation (|82|) it follows that , drPx decay and 



2 6gi|P-| ^ 12gi|P-] 

as w — )■ oo. Hence there is an Z = 2 instability where transverse derivatives of these gauge 
invariant perturbations blow up linearly on the horizon. 

Similar statements for / = 1 even and / > 1 even perturbations can be deduced from table 
[1] and Appendix [Bl The associated functions W in ( I79|) for each variable ^ are easily read 
off from [211 [25] and [29]; in table [1] we provide Vr|r.=M which as shown above is enough to 
determine the existence of a tower of conservation laws and associated instabilities on the 
horizon. 
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A Numerical method 
A.l Algorithm 

Here, we explain our numerical method to solve Eq. fl25|) . We take the domain of the 
numerical calculation as [/ G [—0.5,0] and v G [0,2000]. We discretize coordinates U and 
V as {Ui,Vj) (z = 0,1,2, ■■■ ,iV, j = 0,1,2 ■■•M) where Uq = -0.5, Un = 0, Vq = and 
vm = 2000. Denoting = (j){Ui,Vj), Vij = V{Ui,Vj), 6Ui = f/^+i - Ui and 6vj = Vj+i - Vj, 
we obtain the discretized equation for Eq. fl2Sl) as 

XT T X 

+ - 0,,, L^(V^0),^,/2,,+i/2 + Oi6Uf6v,, 6UM^) , (102) 



where (l^0)j+i/2j+i/2 = {Vi+ij+i(pi+ij+i + Vi+ij(pi+ij + Vij+i(t>ij+i + Vij(f)ij)/4:. Solving above 
equation for we have 



3 X„, XTT X„.3 



+ 0{5Uf5vj,5Ui5v']) , (103) 

where (^0) ■+1/2^+1/2 = (^i+i,j0i+i,j + ^ij+i0ij+i + ^i,i0i,i)/4. Note that (pi^ and 0o,j are 
given by the initial condition fl27|) or fl28|) . Thus, using Eq. fll03p . we can determine the {(pij} 
in the whole domain recursively. 

Now, we consider how we should take the numerical grid, {Ui,Vj). In Fig. [19], we depict 
the potential V{U, v) for fixed v slices. We can see that the potential approaches U = and 
tends to be sharp as v increases. This implies that we need smaller grid size as U approaches 
zero. Thus, we take a multigrid generated by following algorithm: Ui = Ui^i — Cih where 
f/jv = and Q = 1, 4, 16, 64, 256 for < i < A^', A^' < i < 2A^', 2A^' < i < 3A^', 3A^' < i < 4A^' 
and 4A^' < i < 5N', respectively. Here, we define A^' = A^/5 and h = 0.5/ YliLo^ ^i- For this 
choice of h, we have Uq = —0.5. In our calculation, we take the grid number as A^ = 4 x 10^. 
By this choice of numerical grid, we can resolve the potential V{U,v) within v < 2000. For 
v-direction, we take uniform grid 6vj = 0.1. 

A. 2 Evaluating transverse derivatives on the horizon 

In sections [3] and [5l we evaluate 9"0|r=M {n = 1,2,3,4) from the numerical solution 
(f){U, v). The numerical differentiations tend to lose the accuracy as n increases. Thus, we use 
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Figure 19: The potential V{U, v) for m = and / = 0. We take fixed v slices: v = 20, 40, 60. 
We can see that the potential approaches U = and tends to be sharp as v increases. 



a trick to evaluate the differentiations. In the (w, r)-coordinates, the wave equation fl25l) can 
be written as 

F'{r) l{l + i; 



2d^dr(t) + F{r)df.(t) + F'{r)dr(t) - V{r)(t) = , V{r) 

I I - 

where ' = d/dr. From the equation, we obtain recursion equations for 



m 



\r=M 



as 



2d^dr(j)\r=AI 
^2, 



[V(j)]r=M , 

2 



\r=M 



\r=M 



M2 



dr<j) + Vdrd + V'c 



6 



12 

M3 



OA + Vd: 



r=M 
2i , 



2V'dr(j) + V% 



r=M 



2d^d^(j)\r=M 



12 

M2 



48 ^2 



72 



dr 



M4 

+ Vdf(f) + W'dlct) + W'drcj) + V" 



r=M 



(104) 



(n = l,2,---) 

(105) 
(106) 

(107) 



:io8) 



These equations are obtained by operating 5™ and setting r = M in Eq. fll04p where m 
is taken to be m = 0,1,2,3. Once we know 4'\r=M from the numerical calculation, we can 
determine dr(p\r=M by integrating Eq. (11051) with respect to v. From the (p\r=M and dr(p\r=M, 
we can also determine d'^(f)\r=M by integrating Eq. (11061) . In a similar way, we can determine 

dl!^(f)\r=M- 



B Regularity of Moncrief 's variables 

In this section we will explain Moncrief 's parameterisation of metric and gauge field per- 
turbations of the RN solution. The essential point is to expand each type of perturbation in 
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a Regge- Wheeler basis for scalar, vector and tensor harmonics on S*^ which we have written 
covariantly in fl75| 17^ 1771) . 

As explained Moncrief works in the Hamiltonian formalism. Hence the perturbation vari- 
ables consist of the various metric and gauge field components and their conjugate momenta. 
We will need explicit expressions for 6tt^^ , 6E^ . The former can be computed using [25] 

Sn'' = ^ia^'g'^ - 9''9'')[dtS9ki - VJNi - Vi6N,] , (109) 
where Vj is the connection associated to Qij. The latter can be computed from 

6E^ = -^g''6F,, + (^-^ + W - g^'E^Sg^i (110) 

which is valid if the background solution satisfies Ni = 0, Fij = 0, as is the case for the RN 
solution we will consider. Indeed the RN background is given by equations (!73|) and (!74|) . 

We will also show that the subset of Moncrief's variables which appear in the reduced 
gauge invariant perturbation equations, see table [H are smooth on the horizon. 

B.l Odd perturbations 

For odd perturbations Moncrief's parameterisation of the perturbed metric and gauge field 
is [301123] 

Sgijdx'dx^ = 2hidreadx"- + h2eabdx'^dx^ (111) 
SNidx' = hoCadx" (112) 
6Aidx' = Acadx"" . (113) 

The momenta conjugate to 5gra = hiCa and 6gab = h2eab are given by (see e.g. [30] ) 

pi singe" 2p2 singe"'' 

=WTTy' =/(/ + l)(/-l)(/ + 2) ^^'^^ 

where the normalisations are chosen so that pi,p2 is conjugate to hi, /i2 when one integrates 
the action over S'^. The momentum conjugate to 5Aa = Aca is 

Hence perturbations are described by the set of function {A,E, hi, /i2,Pi,P2, ^o)- 

For / > 1, Moncrief performs a canonical transformation to a new set of canonical variables 
{hi, h2, A,pi,p2, E) H- {ki, k2, fi, ni, tt2, tt/), where in particular Hf = A and vti = pi. All the 
new variables are gauge invariant except k2 which is cyclic; its conjugate momentum 712 is 
the only constraint function (which follows by variation with the respect to the lapse ho). 
Hamilton's equations for the remaining pairs (A;i,7ri) and (/i,7rj) can be combined to give 
two second order coupled equations for tt/ and Hg = ttti — 2QI{1 + l)7r//r. These equations 
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can be decoupled into two wave equations of the form f[79l) for ip = P± which are two certain 
constant hnear combinations of 7r/,7rg (these are given in [22] and we have written them in 
the extreme case in equation fl^ and (jHSD)- 

For / = 1 there are no tensor perturbations and so we have /i2 = P2 = 0. Moncrief performs 
a canonical transformation {hi, A,pi, E) h-> (/ci, /i, tti, vt/), where in particular Hf = A and 
TTi = |r^pi — 2QA. Variation of the lapse ho in this case gives the constraint vti = 0, and 
together with the fact that ki is cyclic, this implies tti is a constant. Hamilton's equations for 
the remaining pair (/i, vtj) can be written as a single second order wave equation of the form 
fl78|) for = Pf = 7if — 2Q6a/r where 6a = — 7ri/(6M) is a constant by the previous remarks. 

The above discussion shows that to examine smoothness of P± (/ > 1) and tt/ (/ = 1), 
we need to only examine smoothness of A,pi from the original variables. The latter can be 
computed using fll09p and we find 

Pi = l{l + 1) (dthi - drho + \^ , (116) 

giving an explicit expression in terms of metric components. 

A general odd perturbation of the metric regular on H"*" in EF coordinates is 

dg^ydx^dx'^ = 2{h^dv + hrdr)eadx°' + hcabdx'^dx^ (117) 

where hr,hj},h are smooth on the horizon Ti^. Converting to static coordinates gives 



5g^j_ydx^dx^ = 
Hence comparing to Moncrief gives 



h^dt + { hr + ^ ] dr 



Cadx" + he.abdx^dx'' . (118) 



ho = h^ hi = hr + ^ h2 = h (119) 

F 

and so we deduce that /i0)^2 are regular on the horizon, but hi is not. The momentum 
conjugate to hi f lll6p in EF coordinates is therefore 

Pi = l{l + 1) [d^hr - drK + '^h}j (120) 

which is indeed smooth on the horizon. 

A general odd perturbation of the Maxwell field regular on T-i^ in EF coordinates is 

SF^^dx^dx" = {fvadv + fradr) A Cadx" + fYeabdx" A dx^ , (121) 



where /^a, / are smooth. In static coordinates this becomes 

A eadx"" + fYiabdx" A dx'' . (122) 



SF^^dxf'dx" 



fvadt + { fra + ^] dr 



Comparing to Moncrief gives / = — 1/(/ + 1)A, where we have used de = eV^F = — /(/ + l)Ye. 
Therefore we deduce that A is indeed smooth on "H"*". 
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B.2 Even perturbations 

Even perturbations are more complicated so we will not explain the derivation of the 
variables we need. The metric and gauge field are parameterised as 125] 

Sgijdx'dx^ = F-^H2Ydr^ + 2hidrfadx'' + r\Gfab + Kgab)dx''dx'' (123) 
6Nidx' = HiYdr + hofadx" (124) 

6N = -^F^/^HqY (125) 

6Aidx^ = (ai + a2^r)ydr + a2fadx"- (126) 
6Ao = aoY . (127) 

There are analogous expressions for the conjugate momenta Sir'^^ and 6E^] the only one we 
will actually need is 

5E'' = fisinOY . (128) 
We can compute this explicitly using f lllOp 

/i = -r' I SFtrY^l + Q{Ho-H2 + 2K- l{l + l)G) . (129) 

For / > 1 Moncrief transforms {H21 /ii, G) to the variables (/ci, /c2, /cs, defined by 

k, = K + rFG,r - ^Fh k2 = ^{H2-K)- + (130) 

k^ = G ki = hi. (131) 

Moncrief shows that ki, k2 are gauge invariant. He eventually managed to reduce the pertur- 
bation equations to a pair of coupled second order wave equations of the form ( ITHj) for the 
variableqlj 



_ qWi} ^l)(^+ _2) ^^^^^ ^ ^^^2^^ ^ ^ ^yj^^ ^^32) 

H = F-'^^ where F = /i + g/(/ + 1)^:3 (133) 
Ar 

where A = (Z-l)(/ + 2) + 6M/r-4QV^^- It turns out that F is gauge invariant and therefore 
(Qi, -f?^) are both gauge-invariant. The variables appearing in the decoupled equations in this 
case are certain constant linear combinations R± of (Qi, H). Therefore to examine smoothness 
of the final variables R± we need to examine smoothness only of gi, /C3, /i. 

For / = 1 we must have (7 = 0. Moncrief reduces this system to a single wave equation of 
the form fl78|) for the gauge-invariant variable ip = H wherj^ 

H = f\- where k = FH2 - F^'/^dJrF-^/^K) + K- . (134) 

A r 



^■^In Moncrief's paper [25] the variable Qi is called Q; we have renamed this to avoid confusion with the 
electric charge Q. 

^^Here k2 = Fk2 where k2 is the I = 1 variable Moncrief used [25]: we avoided referring to this to avoid 
confusion with k2 for the I > 1 variable. 
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Hence to examine smoothness of H we need to examine only fi,k2- 

A general even perturbation of the metric regular on the future horizon is 



Sg^ydx^dx" = {hy^dv^ + 2hyrdvdr + hrrdr^)Y + 2{hydv + hrdr)fadx"' + r^{Kgab + Gfab)dx"'dx^ 

(135) 

where all component functions are smooth. Converting to static coordinates gives 



5gfj,ydx^dx" 



hm,dt + 2dtdr ( h„r H — ^ 



2 / , , 2hi,j- hyy 

dr ( hrr + + 



Y 



hydt + { hr + — ] dr 



fadx"" + r^iKgab + Gfab)dx^dx' . (136) 



Comparing to Moncrief 's variables shows K, G are the same and 



ho = hy hi = hr + — 

r 



"'VV 



Hi — hyr + 



H2 — Fhrr ~\~ 2hyj- -\- 



(137) 
(138) 



For / > 1 we deduce that the quantities ki, F'^k2^ k^ given in f ll30p . fll3ip are smooth on the 
horizon, and hence gi also is. For / = 1 we deduce that k2 given in fll34p is smooth on the 
horizon. 

A general even perturbation of the Maxwell field in EF coordinates is 



SF^^dx^'dx" = UYdv Adr + {f^adv + fradr) A fadx" 



which in static coordinates becomes 

SF^^dxf'dx" = UrYdt Adr + 



Therefore we find that fll29p is given by 

/i = -r'U + Qi-Fhr 



fvadt -\- I fra ~l~ 



2K, 



fv 



dr 



A fadx'' 



2K - l{l + l)G) 



(139) 



(140) 



(141) 



and hence this quantity is also smooth on the horizon. 

Putting everything together we deduce: the gauge-invariant functions R± for / > 1 and H 
for / = 1 are indeed smooth on "H"*". 
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